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Abstract: A cadmium tungstate crystal boule enriched in 116 Cd to 82% with mass of 1868 g 
was grown by the low-thermal-gradient Czochralski technique. The isotopic composition of cad- 
mium and the trace contamination of the crystal were estimated by High Resolution Inductively 
Coupled Plasma Mass-Spectrometry. The crystal scintillators produced from the boule were sub- 
jected to characterization that included measurements of transmittance and energy resolution. A 
low background scintillation detector with two 116 CdW04 crystal scintillators (586 g and 589 g) 
was developed. The detector was running over 1727 h deep underground at the Gran Sasso Na- 
tional Laboratories of the INFN (Italy), which allowed to estimate the radioactive contamination of 
the enriched crystal scintillators. The radiopurity of a third 116 CdW04 sample (326 g) was tested 
with the help of ultra-low background high purity germanium 7 detector. Monte Carlo simulations 
of double /3 processes in 116 Cd were used to estimate the sensitivity of an experiment to search for 
double j3 decay of 116 Cd. 
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1. Introduction 

The neutrinoless mode of the double beta decay (0v2j8) is a powerful tool to investigate properties 
of neutrino and weak interactions. Even negative results of double j8 decay experiments provide 
important information on the absolute scale and hierarchy of the Majorana neutrino mass, lepton 
number conservation, right-handed admixtures in weak interaction, existence of majorons, and 
other effects beyond the Standard Model 



- 1 - 



The cadmium 116 is one of the most prospective isotopes to search for 0v2j3 decay thanks 
to the high energy of the decay Qpp = 2809 ±4 keV j^], the relatively large isotopic abundance 
7.49 ±0.18 % [||], the promising theoretical predictions [Q, H, ^] and the existence of cadmium 
tungstate (CdWOzt) crystal scintillators allowing to realize a calorimetric ("source = detector") 
experiment. 

The most sensitive 2/3 experiment to search for double beta decay of 116 Cd was performed 
in the Solotvina Underground Laboratory with the help of cadmium tungstate crystal scintillators 
enriched in 116 Cd ( 116 CdW04, see [f7|] and references therein). The two neutrino mode of 2j8 decay 
of 116 Cd was observed with the half-life Ty 2 = (2.9+°;*) x 10 19 yr, and the half-life limit on Ov2/3 
decay of 116 Cd was set as Tin > 1.7 X 10 23 yr at 90% C.L. It corresponds to an upper bound on the 
effective Majorana neutrino mass (m v ) < 1.7 eV tfjty. Searches for double j3 processes in 106 Cd, 
108 Cd and 114 Cd were realized by using low background CdW04 crystal scintillators produced 
from cadmium of the natural composition [|8[ |9|]. Recently, a cadmium tungstate crystal scintillator 
enriched in 106 Cd was developed JT0[], and an experiment to search for double beta processes in 
106 Cd with the help of this detector is in progress in the Gran Sasso Underground Laboratories 
[11]. In addition, cadmium tungstate crystal scintillators were also successfully applied to study 
the fourth-forbidden j5 decay of 113 Cd ^ [l3|] and to detect, for the first time, the a decay of 180 W 
with the half-life T 1/2 = (l.l±g;|) x 10 18 yr JJ. 

High sensitivity double /3 experiments require detectors with maximal concentration of the 
studied isotope, high detection efficiency to double j8 processes, low (ideally zero) level of ra- 
dioactive contamination and ability of particle discrimination to further reduce the background, 
good energy resolution, large mass, and stability of operation over long (several years) time. 

In the next Section we describe the development of large volume cadmium tungstate crystal 
enriched in 116 Cd. The characterization of 116 CdW04 crystal scintillators produced from the crystal 
is presented in Section 3. The low background detector with the enriched scintillators is described 
in Section 4. The data of low background measurements, both in scintillation mode and with the 
help of ultra-low background high purity germanium (HPGe) y ray spectrometry, were analyzed 
to estimate radioactive contamination of the crystal scintillators. Finally we have simulated the 
detector response to the double j8 decay processes in 116 Cd and estimated the sensitivity to the 
neutrinoless double /3 decay of 116 Cd. 



2. Development of 116 CdWC>4 crystal scintillators 

The production of crystal scintillators from enriched materials should satisfy some specific de- 
mands: minimal loss of expensive isotopically enriched material, high yield of crystals, prevention 
of radioactive contamination. 



2.1 Contamination of enriched 116 Cd 

To produce CdW04 crystals with good scintillation characteristics, it is necessary to control and 
minimize the contamination of the initial materials for the crystal growth at a level of 0.1 — 1 ppm 
for a range of elements. The most dangerous impurities which deteriorate optical and scintillation 
quality of cadmium tungstate crystals are transition metals (Ti, V, Cr, Mn, Fe, Co, Ni, Cu). Con- 
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centration of Ca, Al, Zn, Ag should also be under control. To obtain radiopure scintillators, one 
should avoid contamination by radioactive elements K, Th, U, Ra. 

Samples of enriched cadmium of rather different purity grade were used in the present work. 
Some parts of the enriched material remained after the purification to produce 116 CdW04 crystals 
for the Solotvina experiment, another part of 116 Cd oxide was extracted from the residual after the 
116 CdW04 crystal growth for the experiment tfjfy (see subsection 2.2). One portion of 116 Cd (316 
g) was previously used in experiment to search for double beta decay of 116 Cd to the excited states 



of daughter nuclei [15] 



The contamination of the enriched cadmium was measured with the help of High Resolution 
Inductively Coupled Plasma Mass Spectrometric analysis (Thermo Fisher Scientific ELEMENT2). 
Samples' dissolution was performed by microwave treatment according to the EPA 3052 method. 
Since the solids were not completely dissolved after the microwave digestion treatment, the super- 
natant of samples was analyzed by ICP MS in order to calculate the amount of dissolved sample. 
To enhance the sensitivity to Th and U, the sample solution was treated in an extraction chromato- 
graphic system in order to separate analytes from the matrix. This procedure allowed to reduce 
the dilution factor before analysis up to a value of 150 and to minimize the isobaric interferences 
between Th and WO3 ions. 

Potassium and iron were measured in High Resolution mode. Concentrations were calculated 
based on an external calibration method. We estimate the uncertainties of the measurements as 
about 15% of the given values. Thorium and uranium were measured in Medium Resolution mode 
while other elements were determined in Low Resolution High Sensitivity mode. A semiquantita- 
tive analysis was performed, that is a single standard solution containing some elements at a known 
concentration level (10 ppb of Li, Y, Ce, Tl) was used for calibration and quantification. The uncer- 
tainties when working in semiquantitative mode are about 25% of the given concentration value. 
The contribution of a blank procedure was estimated and subtracted from the samples. The analysis 
results are presented in Table [|. 

2.2 Recovery of enriched cadmium 116 Cd from cadmium tungstate crystalline residue 

A rest after 116 CdW04 crystals growing for the Solotvina experiment was decomposed to extract 
enriched cadmium. Cadmium tungstate is not soluble in acids and alkalis. Molten sodium car- 
bonate (Na2C03, TraceSelect, 99.9999%) was used as a solvent to decompose cadmium tungstate 
crystalline rest. A mixture of 116 CdW04 and Na2C03 in the mass proportion 1 : 1 was prepared, 
which corresponds to approximately 3.3 : 1 molar ratio. Then the mixture was heated in a platinum 
crucible at the temperature 950 °C over 4 hours. As a result the following reaction has occurred: 

Na 2 C0 3 + CdW0 4 -»• Na 2 W0 4 + CdO. 

Sodium tungstate (Na2W04) has the melting point of 696 °C. Therefore this compound was 
in the liquid phase at the temperature of the reaction, while cadmium oxide is insoluble both in the 
molten Na2C03 and Na 2 W04 salts and remained in the form of sediment. After the full decompo- 
sition of the cadmium tungstate the crucible was cooled down to room temperature. Solid Na2C03 
and Na 2 W04 were dissolved in hot deionized water. Finally the precipitation of the enriched cad- 
mium oxide was rinsed and dried. 
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Table 1. Contamination of initial enriched 116 Cd (the average value of several samples) and 116 CdWC»4 
crystal analyzed by ICP-MS analysis. The concentrations of impurities in the samples of 116 Cd oxide and 
the averaged contaminations of initial 116 Cd are normalized on the mass of cadmium. 
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2.3 Purification of enriched cadmium by vacuum distillation and filtering 

Most contaminated samples of enriched 116 Cd were purified by distillation through getters in the 
National Science Center "Kharkiv Institute of Physics and Technology" (Kharkiv, Ukraine) [16]. 



2.4 Synthesis of I16 CdW0 4 compound 

The powder to grow the 116 CdW04 crystal was produced by the NeoChem company (Moscow, 
Russia). All the operations were carried out by using quartz or polypropylene lab-ware, materials 
with low level of radioactive contaminations. Reagents of high purity grade (concentration of any 
metal less than 0.01 ppm) were used. Water, acids and ammonia were additionally distilled by 
laminar evaporation in quartz installation. The high cost of the enriched 116 Cd limits the choice of 
the methods for its additional purification. Recrystallization methods, typically used for the cad- 
mium salts purification, cannot be applied due to the low outcome of the final product (< 85%). 
Therefore, after dissolving the metallic cadmium in nitric acid, the purification was realized by co- 
precipitation on a collector. Additional recrystallization was performed to purify ammonium para- 
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tungstate used as tungsten source. Solutions of cadmium nitrate and ammonium para-tungstate 
were mixed and then heated to precipitate cadmium tungstate: 



Cd(N0 3 ) 2 + (NH 4 ) 2 W0 4 = CdW0 4 + 2NH 4 N0 3 . 

Then the 116 CdW0 4 sediment was rinsed and filtered. Finally the 116 CdW0 4 compound was dried 
and annealed. 



2.5 Growth of 116 CdW0 4 crystal and production of scintillation elements 

The 116 CdW0 4 crystal was grown by the low-thermal-gradient Czochralski technique [jn|, [l9j] 
in a platinum crucible. Under low-thermal-gradient conditions the temperature of the melt was 
close to the melting point of cadmium tungsten crystal, which is « 1270 °C. A crystal boule with 
mass of 1868 g (see Fig. 1, left) was grown from 2139 g of the initial 116 CdW0 4 charge (87% of 
initial charge). 




Figure 1. (Color online) Left: Boule of enriched 116 CdW0 4 crystal. The conic part of the boule is the 
beginning of the crystal growth. Right: Crystal samples cut from the boule: ss 045 x 46.7 mm, 586.2 g, No. 
1 (right); « 045 x 46.1 mm, 589.3 g, No. 2 (middle); « 045.7 x 25.1 mm, 325.6 g, No. 3 (left). 



Three near cylindrical shape crystal scintillators (045 x 46.7 mm, 586.2 g, No. 1; 045 x 
46.1 mm, 589.3 g, No. 2; 045.7 x 25.1 mm, 325.6 g, No. 3) have been cut from the crystal 
boule (see Fig. 1, right). The side surface of the crystals No.l and 2 was diffused with the help of 
grinding paper to reach uniformity of scintillation light collection, which is important to improve 
energy resolution of the detector. 

3. Characterisation of 116 CdWC>4 crystal scintillators 
3.1 Isotopic composition of cadmium in the 116 CdW0 4 crystal 

The isotopic composition of the cadmium in the enriched 116 CdW0 4 crystal was measured with 
the help of the High Resolution Inductively Coupled Plasma Mass Spectrometric analysis. Results 
of the analysis are presented in Table ||[ The absolute isotope abundance for 116 Cd is 82.2%, while 
j3 active 113 Cd has an absolute isotope abundance of 2.14%. The atomic weight of cadmium in the 
116 CdW0 4 crystal is 115.3 as compared to the table value of 112.411 ±0.008 [||]. The atomic 
weight of 116 CdW0 4 molecule is 363.1. 
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Table 2. The absolute isotopic composition of cadmium in the CdWC»4 crystal (%) 



Atomic number Enriched 116 Cd Natural cadmium [pi] 
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12.80 ±0.12 
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28.73 ±0.42 


116 


82.2±0.1 


7.49±0.18 



3.2 Light transmission 

The transmittance of the 116 CdW04 crystal scintillators was measured in the spectral range 330 — 
700 nm using a PERKIN ELMER UV/VIS spectrometer Lambda 18. A thin (2.6 mm) sample of 
the 116 CdW0 4 crystal was placed in the reference beam of the instrument to correct the reflection 
losses. The results of the optical transmission measurements for the 116 CdW04 crystal No. 2 
shown in Fig. |2| demonstrate that the material exhibits reasonable transmission properties in the 
relevant wavelength range 400 — 600 nm of the CdW04 emission spectrum. 




Figure 2. (Color online) The optical transmission curve of CdWC<4 crystal No. 2 before and after anneal- 
ing measured with 2.6 mm sample in reference beam. 

From the data of the transmission measurements we have derived the attenuation length of the 
material 31 ±5 cm at the wavelength of CdW04 emission maximum 480 nm [0]. At 400 nm the 
attenuation lengths of the samples are 5.8 cm (No. 1), 5.2 cm (No. 2) and 4.1 cm (No. 3). One 
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can explain the decrease of transmittance for the samples distant from the beginning of the crystal 
growth by the well known effect of defects increase during CdW04 crystal growth. 

After the low background measurements (see below subsection 4.2) the crystals No. 1 and 2 
were annealed at the temperature m 870 °C over 55 hours. The annealing improved transmittance 
of the samples in the region of wavelengths 350 — 420 nm on 10 — 40 % (see Fig. |2|). 



3.3 Energy resolution 

To measure the scintillation properties, the samples No. 1 and 2 were optically coupled with the 
help of Dow Corning Q2-3067 optical couplant to 3" photomultiplier (PMT) Philips XP2412. To 
improve scintillation light collection, the crystals were wrapped by a few layers of polytetrafluo- 
roethylene (PTFE) tape. The measurements were carried out with 10 /is shaping time of ORTEC 
575 spectroscopy amplifier to collect most of the charge from the anode of the PMT. The detectors 
were irradiated by 7 quanta of 137 Cs, 207 Bi and 232 Th sources. Fig. [3] shows the pulse amplitude 
spectra measured with the 116 CdW04 crystal scintillator No. 2. The energy resolution 11.1% and 
10.1% (FWHM) were obtained for the 662 keV 7 line of 137 Cs with the detectors No. 1 and 2, 
respectively. The energy resolution for 2615 keV 7 line of 232 Th source is 7.1% and 6.7% for the 
scintillators No. 1 and 2, respectively. 
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Figure 3. The energy spectra of 137 Cs, 207 Bi and 232 Th 7 quanta measured by the 116 CdWC>4 scintillation 
crystal No. 2. 



After the annealing at high temperature the energy resolution of the crystals was improved 
(see subsection 4.2). 
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4. Low background measurements, results and discussion 



4.1 Measurements with ultra-low background HPGe 7 ray spectrometry 

The 116 CdW0 4 crystal No. 3 was measured for 788 h with the ultra-low background HPGe 7 ray 
spectrometer GeCris (volume 468 cm 3 , 120% relative efficiency). The background data were accu- 
mulated over 1046 h (see Fig. Q). In order to determine the radioactive contamination of the sam- 
ple, the detection efficiencies were calculated using Monte Carlo simulation based on the GEANT4 



software package [22]. Peaks in the spectra are due to the naturally occurring radionuclides of the 
uranium and thorium chains and 40 K. Only upper limits could be obtained for corresponding activ- 
ities. We have detected low contamination by 137 Cs and 207 Bi in the crystal on the level of 2.1(5) 
mBq/kg and 0.6(2) mBq/kg, respectively. In addition, we have observed peaks of 44 Ti (67.9 keV 
and 78.4 keV, the half-life of 44 Ti is 60 yr) and its daughter 44 Sc (1157.0 keV) in the data accu- 
mulated with the 116 CdW04 crystal. However, the peaks are due to contamination of the HPGe 
detector (not the crystal scintillator sample) by 44 Ti. Indeed the detector, before the run with the 
116 CdW04 sample, has been used to measure a sample of titanium with rather high activity of 
^Ti. Limits on mBq/kg level were obtained for other potential contaminations. The results of the 
measurements are presented in Table B. 
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Figure 4. (Color online) Energy spectra measured with the 325.6 g 116 CdWC>4 sample over 788 h and with- 
out sample over 1046 h (BG) by ultta-low background HPGe 7 spectrometer. Some excess of the specttum 
accumulated with the 116 CdWC>4 sample is due to an accidental contamination of the HPGe detector by 
radioactive 44 Ti. The energy of the 7 lines are in keV. 



4.2 Low background detector system with 116 CdW0 4 crystal scintillators 

The 116 CdW04 crystal scintillators No. 1 and 2 were fixed inside the cavities 047 x 61 mm in 
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the central part of the polystyrene based plastic scintillator light-guides (UPS923A, Amcrys-H, 
Ukraine), 70 mm in diameter and 194 mm in length (a schematic view of the set-up is presented in 
Fig. ||). The cavities were filled with liquid scintillator (LS-221, Institute for Scintillation Materials, 
Kharkiv, Ukraine) which does not affect the polystyrene scintillator. The scintillating light-guides 
act as active veto. A significant difference of CdW04 pulse-shape (effective average decay time is 



13 /is [21]) in comparison to much faster plastic and liquid scintillators signals (few nanoseconds) 
offers the possibility to exploit the discrimination of the light-guide signals. Each plastic light- 
guide was optically connected on opposite sides to two high purity quartz light-guides 070 x 200 
mm each. Two low radioactive 3" diameter PMTs Hamamatsu R6233MOD viewed each detector 
from opposite sides. The light-guides are wrapped by a few layers of PTFE reflection tape. All 
the optical contacts between the light-guides and PMTs were provided by Dow Corning Q2-3067 
optical couplant. The detectors with the 116 CdW04 crystals were placed between two polystyrene 
based plastic scintillators (UPS 89, Amcrys-H, Ukraine) 500 x 300 x 120 mm. A channel 014 x 200 
mm was made (in the middle of the upper plastic scintillator of its width, 5 1 mm from the above 
edge of the plastic) to insert radioactive sources. Two low background 3" diameter PMTs ETL 
9302FLA were optically connected to the plastic scintillators. 




A-A 



Figure 5. (Color online) Low background scintillation detector. (1) 116 CdWC>4 crystals scintillators, (2) 
plastic scintillator light-guides, (3) cavities in the light-guides filled by liquid scintillator, (4) quartz light- 
guides, (5, 6) plastic scintillators, (7) four PMTs Hamamatsu R6233MOD, (8) two PMTs ETL 9302FLA, 
(9) channel 014 x 200 mm in the plastic scintillator 6 to insert radioactive sources. 

The detector system was installed deep underground in the low background DAMA/R&D 
set-up at the Gran Sasso National Laboratories of the INFN (Italy). The detector system was 
surrounded by Cu bricks and sealed in a low radioactive air-tight Cu box continuously flushed 
by high purity nitrogen gas (stored deeply underground for a long time) to avoid the presence 
of residual environmental radon. The Cu box was surrounded by a passive shield made of high 
purity Cu, 10 cm of thickness, 15 cm of low radioactive lead, 1.5 mm of cadmium and 4 to 10 
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cm of polyethylene/paraffin to reduce the external background. The shield was contained inside a 
Plexiglas box, also continuously flushed by high purity nitrogen gas. 

An event-by-event data acquisition system based on a 1 GS/s 8 bit transient digitizer (Acqiris 
DC270) records the time of each event and the pulse shape over a time window of 100 /^s from 
the 116 CdW0 4 detectors (the sum of the signals from two PMTs), the plastic scintillator, and the 
sum of signals from the 116 CdW04 scintillators attenuated to provide an energy scale up to « 10 
MeV (the electronic chain of the detector system is briefly summarized in Fig. ^). Taking into 
account the slow kinetics of the CdW04 scintillation decay, the sampling frequency of the transient 
digitizer was set to 20 MS/s. An especially developed electronic unit (SST-09) provides triggers for 
cadmium tugstate scintillation signals. The unit rejects PMT noise, plastic scintillator light-guide 
signals and CdW04 events with large admixture of the plastic. Further rejection of the plastic 
overlaps can be realized off line by the pulse-shape analysis described below. 
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Figure 6. (Color online) Schema of the electronic chain (details of the low background detector are de- 
noted in Fig. 5). (FAN IN/OUT) Linear FAN-IN/FAN-OUT, LeCroy Model 428F; (SST-09) home made 
electronic unit to provide triggers for cadmium tugstate scintillation signals; (Attenuator) Dual attenuator, 
CAEN model No 110; (PC) personal computer. Inputs and outputs of the electronic units are denoted as "In" 
and "Out", respectively; "Trg" denotes a trigger input of the transient digitizer. 

The energy scale and resolution of the detector system was tested with 22 Na, 60 Co, I37 Cs and 
228 Th 7 sources. The energy resolution of the 116 CdW04 detectors before the annealing can be de- 
scribed by functions: FWHM r i = ^10.5 x E y , and FWHM^ = ^540 + 8.5 x E y for the detectors 
No. 1 and 2, respectively (here E y and FWHM r are in keV). The energy scale and resolution of 
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the detectors were tested once during the measurements and at the end of the experiment with the 
help of 22 Na 7 source inserted into the set-up through the special channel without switching off the 
high voltage of the PMTs. We have observed neither shift of the energy scale nor degradation of 
the energy resolution of the detectors during almost 2.5 months of low background measurements. 

As a result of the annealing of the 116 CdW04 crystal scintillators performed after the low 
background measurements, the energy resolution of the detectors for the 2615 keV / line of 228 Th 
was improved from 6.9% to 5.3% for the detector with the crystal No. 1, and from 6.2% to 5.0% 
for the detector No. 2. Energy spectra accumulated with the 116 CdW04 detector No. 2 after the 
annealing are presented in Fig. [7|. 




Figure 7. (Color online) Energy spectra accumulated by the 116 CdWC<4 detector No. 2 with 137 Cs, 60 Co 
(a), 22 Na and 228 Th y sources (b) in the low background set-up after annealing of the crystals. Energies of y 
lines are in keV. 



4.3 Data analysis 

The energy spectra accumulated with the 116 CdW04 detectors over 1727 h of low background 
measurements are presented in Fig. |8[ The spectra are rather similar with a small difference in the 
region 0.6 — 1.2 MeV where a particles from U/Th are observed (see Subsection 4.3.1 describing 
a pulse-shape discrimination to select a particles). The counting rate of 0.26 count/s in the energy 
interval 0.08-0.6 MeV is mainly due to the decay of 113 Cd (Q p = 320 keV, T l/2 = 8.04 x 10 15 yr) 
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with the activity (0.10±0.01) Bq/kg 1 and 113m Cd (Qp = 583 keV, T l/2 = 14.1 yr) with the activity 
(0.46 ±0.02) Bq/kg. 
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Figure 8. (Color online) The energy spectra accumulated with the 116 CdWC<4 crystal sintillators No. 1 
and 2 in the low background DAMA/R&D set-up over 1727 h. The energy spectra of a events selected by 
the pulse-shape discrimination (see text) are also shown. Some difference in the data in the energy region 
0.6—1.2 MeV visible both in the raw data and in the spectra of a events is due to higher a activity of U/Th 
daughters' traces in the crystal No. 2. In the inset, the a spectrum of the detector No. 2 is depicted together 
with the model, which includes a decays from 238 U and 232 Th families. 

Contributions to the background above the energy 0.6 MeV were analyzed by the pulse-shape 
discrimination and the time-amplitude techniques. 

4.3.1 Pulse-shape discrimination 

To select / quanta (j8 particles) and a particles, the data of the low background measurements were 



analyzed by using the optimal filter method proposed by E. Gatti and F. De Martini [23] (see also 
j p4[ ] where the analysis was developed for CdW04 crystal scintillators). For each experimental 
signal, the numerical characteristic of its shape (shape indicator, SI) was defined as 5/ = x 
P{h) /Hf{tk), where the sum is over time channels k, starting from the origin of the signal and up 
to 30 /is, f(tk) is the digitized amplitude (at the time tt) of a given signal. The weight function P(t) 
is defined as: P(t) = {f a {t) — fy(t)}/{f a (t) +fM)}, where the reference pulse shapes f a (t) and 
fy{t) are taken from The pulse-shape discrimination of the events accumulated in the low 
background measurements with the 116 CdW04 detector No. 1 is demonstrated in Fig. ^. Alpha 
events were selected from the accumulated data (see the energy spectra of a particles in Inset of 



'The activity is calculated on the basis of the isotopic abundance and the half-life of I13 Cd. 
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Fig. ||), which allow to estimate the total internal a activity of U/Th as 1.9(2) mBq/kg and 2.7(3) 
mBq/kg in the 116 CdW04 crystal scintillators No. 1 and 2, respectively. Slightly higher a activity 
in the sample No. 2 can be explained by accumulation of U/Th daughters' traces in the melt during 
the crystal growth. 

To estimate the activity of the a active nuclides from the U/Th families in the crystals, the 
energy spectra of the a events were fitted in the energy interval 0.5 — 1.5 MeV by using a simple 
model: ten Gaussian functions to describe a peaks of 232 Th (and its daughters: 228 Th, 224 Ra, 
220 Rn, 216 Po, 212 Bi), 238 U (and its daughters: 234 U, 230 Th, 210 Po) plus exponential function (to 
describe background). Fit of the a spectrum accumulated with the detector No. 2 is shown in Inset 
of Fig. ||. Because of the worse energy resolution for a particles in comparison to y quanta in 
CdW04 scintillation detectors [14], we conservatively give limits on activities of 232 Th and 238 U, 
230 Th, 210 Po (expected to be not in equilibrium with 238 U) in the 116 CdW04 scintillators. The data 
obtained from the fits are presented in Table |5[ 

In addition, the optimal filter method allows to distinguish and reject from the data also over- 
lapping of plastic scintillator's pulses with 116 CdW04 signals, random coincidence of events, some 
events from the fast chain of decays 212 Bi— s> 212 Po— > 208 Pb ( 232 Th family). 
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Figure 9. (Color online) Shape indicators (see text) versus energy for background exposition over 1727 h 
with 116 CdW04 crystal scintllator No. 1 in the low background set-up. Three sigma intervals for shape 
indicator values corresponding to y quanta (j3 particles) and a particles are drawn. 

4.3.2 Time-amplitude analysis of 228 Th and 227 Ac contaminations 

Activities of 228 Th ( 232 Th family) and 227 Ac ( 235 U) in the 116 CdW0 4 crystal scintillators were 
estimated by the time-amplitude analysis of the events accumulated in the low background mea- 
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surements. The technique of the time-amplitude analysis is described in detail in [0, |25], [26| ]. 



To determine the activity of 228 Th, the following sequence of a decays was selected: 224 Ra 



224t 



{Qa = 5.789 MeV, T 1/2 = 3.66 d) -> 220 Rn (Q a = 6.405 MeV, T 1/2 = 55.6 s) 21b Po (Q a = 6.907 



216t 



MeV, Tin = 0.145 s) — > 212 Pb. The obtained a peaks and the distributions of the time intervals 



between events (see Fig. |10() are in agreement with those expected for a particles of the chain. 
Taking into account the efficiencies in the time windows to select 22() Rn — } 216 Po — > 212 Pb events 
(94.5%), the activities of 228 Th in the crystals No. 1 and 2 were calculated as 0.057(7) mBq/kg and 
0.062(6) mBq/kg, respectively. 

By using positions of the three a peaks in the y scale of the detector we have obtained the 
following dependence of a/fi ratio on energy of a particles E a : a/fi = 0.113(6) +0.132(10) x 
10~ 4 £ a in the energy interval 5.8 — 6.9 MeV (E a is in keV). 

The same approach was used to search for the chain 223 Ra (Q a = 5.979 MeV, 7\/ 2 = 11.44 d) 
219 Rn (Q a = 6.946 MeV, T l/2 = 3.96 s) 215 Po (Q a = 7.526 MeV, T l/2 = 1.78 ms) -> 211 Pb 
from the 235 U family. Limit on activities of 227 Ac in the crystals No. 1 and 2 was obtained on the 
level of < 0.002 mBq/kg. 
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Figure 10. (Color online) Alpha peaks of 224 Ra, 220 Rn and 216 Po selected by the time-amplitude analysis 
from the data accumulated during 1727 h with the 116 CdWC<4 detector No. 1. The obtained half-lives 
of 220 Rn (41+^ s) and 216 Po (0.162+JJ;^ s) are in agreement with the table values (55.6 s and 0.145 s, 
respectively P7P). 



4.3.3 Selection of Bi-Po events 

Double pulses from the decays of 214 Bi {Qp = 3.272 MeV) -> 214 Po (Q a = 7.833 MeV, T l/2 = 164 
— > 210 Pb (daughters of 226 Ra from the 238 U family) were searched for. One event in the data 
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accumulated with both detectors was found. Therefore, taking into account the efficiency of the 
Bi-Po events selection in the time interval of the analysis 1 — 75 /is (26.7%), we set limit on 226 Ra 
activity in the crystals No. 1 and 2 as < 0.005 mBq/kg. 

A summary of the radioactive contamination of the 116 CdW04 crystal scintillators is presented 
in Table H 

Table 3. Radioactive contamination of 116 CdW0 4 crystal scintillators. Data for CdW0 4 @, 116 CdW0 4 
[R fj\, and 106 CdWC>4 O] are given for comparison. 



Chain Nuclide 



Activity (mBq/kg) 





(Sub-chain) 




116 CdW0 4 




CdW0 4 


116 CdW0 4 


106 CdWO 4 






INO. 1 


INO. Z 


INO. J 
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iff21 rTh 
|/p 






Jv 






\ zo 


<-- c 


U.-5(l) 


\ 1 1 




CO 






\ U.4 / 


\ U.4 








110m a „ 
A g 


U.Uu(4; 


a c\f.i a \ 

U.Uo(4 ; 












Cd 


1 AA/ 1 A\ 

1UU(1U) 


1 AA / 1 A\ 

1UU(1U) 




558(4) 


91(5) 


1 HA 
1 /4 




113m Cd 


460(20) 


460(20) 




< 3.4 


1.1(1) 


112000(5000) 




137 Cs 






2.1(5) 


<0.3 


0.43(6) 






207 Bi 






0.6(2) 






1.3(3) 


232 Th 


232 Th 


<0.08 


<0.08 




< 0.026 


0.053(9) 


<0.1 




228 Th 


0.057(7) 


0.062(6) 


< 2.0 


0.008(4) 


0.039(2) 


0.053(5) 


235TJ 


235tj 






<4.0 










227 Ac 


< 0.002 


< 0.002 




0.014(9) 


0.0014(9) 




238TJ 


238TJ 
234m p a 


<0.4 


<0.6 


<58 


< 0.045 


<0.6 
<0.2 


<0.3 




230 Th 


<0.06 


<0.05 






<0.5 


<0.8 




226 Ra 


< 0.005 


< 0.005 


< 2.6 


< 0.018 


< 0.004 


<0.3 




210p b 






< 15000 




<0.4 






210p o 


<0.4 


<0.6 




< 0.063 







Total a activity 1.9(2) 2.7(3) 



0.26(4) 1.40(10) 



<0.3 
2.1(1) 



Selection of double pulses produced by the fast chain of the decays 212 Bi (Qp = 2.254 MeV) 
-> 212 Po (Q a = 8.954 MeV, Ty 2 = 0.299 /is) -> 208 Pb ( 232 Th family) was developed. An example 
of such an event pulse-shape and the result of the selection for the 116 CdW0 4 crystal scintillator 
No. 1 are presented in Fig. 11, The analysis gives the activity of 212 Bi (which is in equilibrium 
with 228 Th) 0.054(5) mBq/kg (crystal No. 1) and 0.095(6) mBq/kg (crystal No. 2) in a reasonable 
agreement with the results of the time-amplitude analysis. 
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Figure 11. (Color online) (a) Example of 212 Bi^ 212 Po^ 208 Pb event in the 116 CdW0 4 scintillator, (b) The 
time distribution for the fast sequence of /3 ( 212 Bi) and a ( 212 Po) decays selected by the pulse-shape and the 
front edge analyses from the background data accumulated with the 116 CdWC<4 detector No. 1 over 1727 h. 
The fit of the time distribution gives a half-life 7jy 2 = (0.29 ±0.03) jxs which is in good agreement with the 



table value for 212 Po (0.299 fis Wf\). 



4.4 Background in the region of Qpp of 116 Cd 

By using the pulse-shape and the front edge analyses we can substantially reduce the background 
of the detector in the energy region of interest near 2.8 MeV where a peak of neutrinoless double /3 
decay of 116 Cd is expected (see Fig. 12). The counting rate of the detector after removing the part of 



data accumulated during first 17 days of measurements (to avoid effects of radon and cosmogenic 
activation) in the energy interval of interest 2700 — 2900 keV is 0.28 counts/(yrxkeVxkg). 

A certain contribution to the background « 0.1 counts/(yrxkeVxkg) in the energy interval 
2.7 — 2.9 MeV comes from the 2615 keV peak caused by external /quanta from the decays of 208 T1 
(daughter of 232 Th). Contamination of the set-up (first of all of PMTs, cables, quartz light-guides, 
copper shield) by thorium can be a source of the background. We are going to simulate 208 T1 decays 
to estimate the contribution from different parts of the set-up. Our preliminary calculations show 
that the PMTs and quartz light-guides can be the main sources of the 2615 keV peak. We are going 
to apply ultra-low background HPGe y ray spectrometry to measure radioactive contamination of 
the PMTs and sample of the quartz. 

The two neutrino mode of 2/3 decay of 116 Cd (assuming a half-life T?J 2 = 2.8 x 10 19 yr [|28|]) 



1/2 

contributes 0.00045 counts/(yrxkeVxkg) in the energy interval 2.7 — 2.9 MeV. 

One of the major sources of the detector background in the energy region of the expected Ov2/3 
peak of 116 Cd is contamination of the 116 CdW04 crystals by 208 T1. The Monte Carlo simulation 
of internal 208 T1 (see Fig. [l2|, c) gives about sa 0.09 counts/(yrxkeVxkg) in the region of interest. 
We consider as a possibility the crystals recrystallization to reduce the thorium contamination. 
To our knowledge there is no data in literature about segregation of thorium to CdW04 crystals. 
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Figure 12. (Color online) (a) Initial sum spectrum of the two CdWC»4 detectors measured over 1322 h 
(total exposure 1554 kgxh) in anti-coincidence with the plastic scintillation counter and the active light- 
guides; (b) the spectra of a and j3 + a events selected by using the pulse-shape and the front edge analyses 
(see text) together with the simulated response function for the 212 Bi— ^ 212 Po— > 208 Pb decay chain; (c) /3 and 
7 events selected with the help of the pulse-shape and the front edge analyses (the efficiency of the selection 
procedure for / quanta / j3 particles is 95%). The fit of the 208 T1 7 peak with the energy 2615 keV is shown 
by solid line. The Monte Carlo simulated energy spectra of internal 110m Ag and 208 T1 in the 116 CdW04 
crystals are presented. 



A positive result of recrystallization procedure can be reached if the segregation coefficient for 
thorium is much less than 1. To determine the segregation of thorium in CdWO/t, we are going to 
measure the activity of 228 Th in the 116 CdW04 crystal No. 3 by the scintillation method. Besides 
we are going to measure the activity of 228 Th in the rest of the melt after the crystal growth by 
ultra-low background HPGe / ray spectrometry, and concentration of 232 Th by the High Resolution 
Inductively Coupled Plasma Mass Spectrometry. Higher concentration of thorium in the rest will 
be an indication of low segregation of thorium in CdW04 crystal. Then the crystal No. 3 can be 
recrystallized and measured by the low counting scintillation method to estimate activity of 228 Th. 

In case of positive result all the crystals can be recrystallized. In case of not enough high 
efficiency of the recrystallization to reduce thorium, one could recover the enriched cadmium from 
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the crystals and scraps, purify by physical (vacuum distillation and filtering) 2 and chemical (by 
coprecipitation on a collector) methods, and grow crystal again. Output of purified materials can be 
decreased to reach deeper purification of enriched cadmium. It should be stressed that a satisfactory 
radiopurity level (activity of 228 Th less than 0.01 mBq/kg) was reached in some samples of CdW04 
crystals [M O]. Besides, very high level of radiopurity (~ 0.002 mBq/kg of 228 Th) was detected 



in zinc tungstate (ZnW04) crystal scintillators [29], which have chemical and physical properties 



rather similar to CdW04. An encouraging result of recrystallization was reported recently for 



calcium tungstate (CaW04) crystal scintillators [30]. One could also expect improvement of the 



crystal scintillators quality thanks to deep purification of initial materials and recrystallization. 
Therefore contribution from the 2615 keV peak could be suppressed further. 

It should be also mentioned a "natural" way of the background decrease during next few years 
due to the decay of the trace 228 Th in the 116 CdW04 crystals (assuming broken equilibrium of 232 Th 
chain and lower activity of 228 Ra in comparison to 228 Th and 232 Th). Such an effect (decrease of 
the 228 Th activity in the 116 CdW04 crystals) was observed in the Solotvina experiment [J7|]. 

Finally, we assume that the essential part of the background beyond the 2615 keV peak is due 
to cosmogenic activation of the 116 CdW0 4 crystals by 110m Ag (Qp = 3.0 MeV; T l/2 = 250 d) [0], 
which can provide background up to 3 MeV. Moreover, radioactive isotope 1 10m Ag with activity 



0.4 mBq/kg was observed in [15] in the enriched cadmium 116 used to produce the scintillators. 
Our assumption was justified by the Monte Carlo simulation of the 116 CdW04 detector response 
to internal 11() '"Ag. A simulated model of 110m Ag corresponding to activity 0.06 mBq/kg is shown 
in Fig. 12, c. Another possible cosmogenic radionuclide in the I16 CdW04 crystals can be 106 Ru 



(Qp = 40 keV; T l/2 = 374 d) 106 Rh (Q p = 3.5 MeV; T l/2 = 30 s) 0]. However, cosmogenic 
background is expected to be reduced substantially due to decay of cosmogenic radionuclides, in 
particular of 110m Ag. 

Now the experiment is in progress with an improved energy resolution after the annealing of 
the crystals. We expect to improve the background reduction at this stage of experiment thanks 
to more careful pulse-shape analysis to reject events caused by 228 Th daughters. In particular, we 
use now a higher resolution of the transient digitizer (50 Ms/s instead of 20 Ms/s) to suppress 
further the background caused by the fast sequence 212 Bi — > 212 Po — > 208 Pb. In addition, we are 
developing the analysis of the sequence of the decays 212 Bi (Q a = 6.207 MeV, Tin = 60.55 m) 
208 T1 (Qp = 5.001 MeV, T l/2 = 3.053 m)-> 208 Pb to reject events of 208 T1 decay. Contribution 
from the 2615 keV y line is expected to be reduced to « 0.06 counts/(yrxkeVxkg) thanks to the 
improvement of the energy resolution. 

5. Monte Carlo simulation of double /3 decay of 116 Cd, estimations of experimental 
sensitivity 

The computer simulation of the different radioactive processes in the scintillation low background 
detector with 116 CdW04 scintillators in the 4 71 active shielding has been performed with EGS4 



package [32]. The initial kinematics of the particles emitted in the decay of the nuclei was given by 



2 Unfortunately, only part of the enriched cadmium 116 was purified by the vacuum distillation and filtering (see 
subsection 2.3). 
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an event generator DECAYO [33]. The following double j3 processes in 116 Cd have been simulated: 
0v2j3 and 2v2j3 decay to the ground state and to the five lowest excited levels of 116 Sn; neutrinoless 



double j3 decay with emission of one, two and bulk [34] majoron(s). Approximately 3 — 9 millions 
of decays were simulated for the different channels of 116 Cd 2/3 decay. The calculated distributions 
are shown in Fig. [l3|. In particular, the calculations give the detection efficiency in a peak of 
the neutrinoless double j8 decay of 116 Cd as 89% (one can compare with the value of 83% in the 
Solotvina experiment [Q] where the crystals of the smaller volume were used). 
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Figure 13. (Color online) Simulated response functions of the CdWC<4 detector for the different modes 
of double j3 decay of 116 Cd to the ground state of 116 Sn (a); Simulated distributions of the detector for two 
neutrino (b) and neutrinoless (c) decay of 116 Cd to the excited levels of 116 Sn. Approximately 3 — 9 millions 
of decays were simulated for the different decay channels. The distributions are normalized to 1000 decays, 
except the 0v2j3 mode, which is normalized to 100 decays. 

On the basis of the detection efficiency and of the number of 1 16 Cd nuclei in the two 1 16 CdW04 
crystal scintillators No. 1 and 2 (1.6 x 10 ), assuming decrease of the detector background to 
the level of 0.01 —0.1 counts/(yrxkeVxkg) (expected due to the decay of cosmogenic radionu- 
clides, improvement of the energy resolution and the pulse-shape analysis, possible change of the 
most contaminated parts of the set-up and reduction of 228 Th activity after recrystallization of the 
116 CdW04 crystals), one can estimate the sensitivity of a 5 years experiment to the 0v2j3 decay 
of 116 Cd as Tin ~ (0.5 — 1.5) x 10 24 yr. According to the recent calculations of matrix elements 
[Q, HI, §1, these half-lives correspond to the effective neutrino mass (m v ) 0.4 — 1.4 eV, which is 
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on the level of the most sensitive 2/3 experiments. 
6. Conclusions 

Cadmium tungstate crystal scintillators were developed from cadmium enriched in 116 Cd for a 
high sensitivity experiment to search for double beta decay of 116 Cd. Samples of enriched cad- 
mium were purified by chemical methods, the most polluted part was additionally purified by 
vacuum distillation. Some part of the enriched material was recovered from scraps after enriched 
116 CdW04 crystal growth in 1988. Cadmium tungstate compounds were synthesized from solu- 
tions. A 116 CdW04 crystal boule with mass of 1868 g (which is 87% of the initial charge) was 
grown by the low-thermal-gradient Czochralski technique. Two 116 CdW04 crystal scintillators 
(586 g and 589 g) produced from the boule show an energy resolution FWHM pa 7% (for 2615 keV 
y line of 2()8 T1) in measurements with the crystals directly viewed by photomultiplier. 

The absolute isotopic composition of 116 Cd in the crystal was determined as 82% by mass- 
spectrometry. Thanks to the careful purification of the initial materials and using the low-thermal- 
gradient Czochralski technique, the crystal has rather high transmittance: the attenuation length is 
31 ±5 cm at the wavelength of CdW04 emission maximum 480 nm. After the low background 
measurements the crystals No. 1 and 2 were annealed at the temperature pa 870 °C over 55 hours. 
The annealing improved transmittance of the samples in the region of wavelengths 350 — 420 nm 
on 10-40 %. 

The low background detector system with two enriched 116 CdW04 crystal scintillators (586 
g and 589 g) was installed in the underground Gran Sasso National Laboratories of the INFN 
(Italy). The energy resolution of the detector with 116 CdW04 crystal scintillators inside the plastic 
scintillator light-guides with length of 28 cm was even slightly better than that in the case when 
the scintillators were directly viewed by a photomultiplier. Furthermore the energy resolution was 
improved to Pa 5% (for 2615 keV y line of 208 T1) after the annealing of the 116 CdW0 4 crystals. 

The low background measurements over 1727 h allowed to estimate radioactive contamina- 
tion of the 116 CdW04 scintillators. In addition the radioactive contamination of another 116 CdW04 
sample (326 g) was tested with the help of ultra-low background HPGe 7 ray spectrometry. The ac- 
tivities of 22 Ra and 228 Th, which are the most dangerous isotopes for double j3 decay experiments, 
are on the level < 0.005 mBq/kg and pa 0.06 mBq/kg, respectively. 

By using the pulse-shape discrimination of the 212 Bi— 212 Po events and the anti-coinciden- 
ce signals in the plastic scintillator light-guide and active shield, we have obtained a background 
counting rate of 0.28 counts/(yrxkeVxkg) in the region of interest 2700 - 2900 keV. 

We have tried to estimate the main sources of the detector background on the basis of the 
low background measurements and the Monte Carlo simulation: they are cosmogenic activation 
(most probably 110m Ag), contamination of the 1 CdW04 crystals by 228 Th ( 232 Th family), con- 
tribution of the 2615 keV 7 peak of 208 T1 from details of the set-up. We are going to simulate 
all possible components of background; check radioactive contamination of PMTs, quartz, cables, 
copper (with aim to change the most contaminated elements); recrystallize the crystals to reduce 
the concentration of thorium. 

The low background measurements to search for double j3 decay of 1 Cd with the help of 
the enriched cadmium tungstate crystal scintillators are in progress. The decrease of thorium con- 
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centration in CdW04 crystal scintillators, the decay of cosmogenic nuclides, the improvement 
of the pulse-shape analysis could further decrease the background to the level of ~ 0.01 —0.1 
counts/(yrxkeVxkg), and therefore, improve sensitivity of the experiment up to 7i/ 2 ~ (0.5 — 
1.5) x 10 24 yr over 5 years of measurements. It corresponds, taking into account the recent calcu- 
lations of matrix elements [Q |5| to the effective neutrino mass (m v ) « 0.4 — 1.4 eV. 

Further progress in the experiment to search for the double beta decay of 116 Cd can be ad- 
vanced by applying a massive array of 116 CdW04 crystals with improved energy resolution and 
good particle discrimination ability. 
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